Constraints from the CMB temperature and other common observational data-sets on 

variable dark energy density models 



Philippe Jetzer^* and Crescenzo Tortora^t 
^ Universitiit Zurich, Institut fur Theoretische Physik, 
Winterthurerstrasse 190, CH-8057, Zurich, Switzerland 



o 

(N 



(N 

o 
u 

6 



> 

o 
o 



% 



The thermodynamic and dynamical properties of a variable dark energy model with density scal- 
ing as px oc (1 + z)"*, z being the redshift, are discussed following the outline of Jetzer et al. [l|]. 
This kind of models are proven to lead to the creation/disruption of matter and radiation, which 
affect the cosmic evolution of both matter and radiation components in the Universe. In particular, 
we have concentrated on the temperature-redshift relation of radiation, which has been constrained 
using a very recent collection of cosmic microwave background (CMB) temperature measurements 
up to 2 ~ 3. For the first time, we have combined this observational probe with a set of independent 
measurements (Supernovae la distance moduli, CMB anisotropy, large-scale structure and obser- 
vational data for the Hubble parameter), which are commonly adopted to constrain dark energy 
models. We find that, within the uncertainties, the model is indistinguishable from a cosmological 
constant which does not exchange any particles with other components. Anyway, while temperature 
measurements and Supernovae la tend to predict slightly decaying models, the contrary happens if 
CMB data are included. Future observations, in particular measurements of CMB temperature at 
large redshift, will allow to give firmer bounds on the effective equation of state parameter w^a of 
this kind of dark energy models. 
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I. INTRODUCTION 

The current standard cosmology model relies on the 
existence of two unknown dark components, the so called 
"dark matter" (DM) and "dark energy" (DE), which 
amount to ^ 25% and ^ 70% of the total energy bud- 
get in the Universe, respectively. According to several 
observations, the Universe is spatially flat and in an ac- 
celerated phase of its expansion 0-la|- DE, described as a 
cosmological constant A in its simplest form, is modelled 
by a fluid with a negative pressure, which is a fundamen- 
tal ingredient to explain the actual accelerated expansion 
of the Universe. 

Several models have been proposed to explain DE 
[7H14J. An alternative consists to consider a phe- 
nomenological variable DE de nsity with continuous cre- 
ation/disruption of photons [l|, [fa - KLSJ or matter [l^, 120] ■ 
The DE might decay/grow slowly in the course of the cos- 
mic evolution and thus provide the source/sink term for 
matter and radiation. Different such models have been 
discussed and strong constraints come from very accu- 
rate measurements of the cosmic microwave background 
(CMB) radiation and other typical cosmological probe. 

CMB radiation is the best evidence for an expand- 
ing Universe starting from an initial high density state. 
Within the Friedmann-Robcrtson- Walker (FRW) mod- 
els of the Universe the radiation, after decoupling, ex- 
pands adiabatically and scales as (1 -t- z), z being the 



redshift |2l| . If we assume that each component is not 
conserved, contrarily to the standard scenario, then de- 
pending on the decay mechanism of the DE, the created 
photons could lead to distortions in the Planck spectrum 
of the CMB, and change the evolution of its tempera- 
ture. The chance to appreciate the deviation from the 
standard temperature evolution is given by the increas- 
ingly number of recent works collecting observations of 
CMB temperature both at low [22, 123] and higher red- 
shifts in. 



Following the theoretical lines of |15l - ll8| . in Jetzer ct 
al. yj we have discussed a variable DE model A(z) ex 
(1 4- z)™ decaying into photons and DM particles. In 
particular we studied thermodynamical aspects in the 
case of a continuous photon creation, which implies a 
modified temperature redshift relation for the CMB. We 
have tested the predicted temperature evolution of radi- 
ation with some recent data on the CMB at higher red- 
shift from both Sunyaev-Zcl'dovich (SZ) effect and high- 
redshift QSO absorption lines. 

In this paper we will present the results obtained in [l| 
within a more detailed theoretical background, by fur- 
ther discussing the main hypothesis which concern the 
energy conservation and thermodynamical laws. Further- 
more, wc will test our model with an updated collection 
of data, joining our sample with the new temperature 



observations reported in |24| . For the first time we will 
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combine to the temperature measurements also the data 
from different kind of observations, like distance moduli 
of Supernovae la, observations of the CMB anisotropy 
and the large-scale structure, together with observational 
Hubble parameter estimations. We will dedicate a partic- 
ular attention to the estimate of both the actual matter 



density parameter, the DE parameter m and the effective 
equation of state parameter Weg. The impact of each ob- 
servational probe on these estimates is also investigated. 
The paper is organised as follows. In f|ll]wc will present 
the model, while the data sets and fitting procedure are 
described in mill The results are discussed in ijIVI and JfVl 
is devoted to a discussion of the results and conclusions. 



II. THEORETICAL FRAMEWORK 



A. Friedmann equations 



with w'' being the four velocity. After easy calculations, 
this conservation law reads 



Ptot + Hptot +Ptot)H = 



(2.4) 



where H = R/R is the Hubble parameter. In the stan- 
dard approach, each component is conserved, thus Eq. 
(|2.4p holds for all the fluid components, but here we will 
suppose that each component will exchange energy with 
each other. In particular, for matter, radiation and DE 
we impose the following relations (see also |25| ) 



Pri 



3p,nH = (1 - e) C, 



(2.5) 



We assume a cosmological framework based on the 
usual Robertson- Walker (RW) metric element [2l[ and 
that the Universe contains three different components: 
a) a matter (both baryons and DM) fluid, with equation 
of state pm = (since Pm << Pm), b) a generalised fluid 
with pressure p-y = (7— l)p-y, where 7 is a free parameter, 
which is set to 4/3 for a properly said radiation fluid, and 
c) a DE, X component, with pressure p^ and density px- 
The equation of state for the x component could assume 
a very general expression, but we limit ourselves to con- 
sider the simple linear relation p^ = WxPx- We will set 
any 'bare' cosmological constant Ag equal to [l|. With 
these components we get for the Einstein field equations 
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where ptot ^ Pj+Px and ptot ^ Pm + Pj + Px are the total 
pressure and density, R is the scale factor, fc = 0, ±1 is 
the curvature parameter and a dot means time derivative. 
Furthermore, we will assume that there is no curvature, 
thus fc = 13. 

In the following we will adopt Wx = —1, however 
since we are assuming that the vacuum decays into radia- 
tion and massive particles, the effective equation of state 
Weff , which is the measured quantity, can differ from — 1 
(see later for further details). 



B. Cosmological conservation law^s and density 
evolution 

Following Jetzer ct al. [1[ (but see also [l^ [T^l) we 
set here the energy conservation equation for the differ- 
ent fiuids. The fiuid as whole, verifies the Bianchi iden- 
tity, which means that energy and momentum are locally 
conserved. In formulae it is V^T'"' = 0, with T^"' the 
stress-energy tensor 



Pj + S^PyH = e Cx , 
Px + i{Px + Px)H = -Cx 



(2.6) 
(2.7) 



where we assume that both the matter and the radiation 
fluids exchange energy with the DE as parameterized by 
Cx and e. In particular, Cx depends on the DE and, 
indeed, acts as a source/sink term for the fluids energy. 
Evidently, if no interaction between the different com- 
ponents exists, then Cx is null and the standard picture 
is recovered. Moreover, if e := (e = 1), then the DE 
exchanges all the energy with matter (radiation). Cx 
can describe different physical situations such as, for in- 
stance, a thermogravitational quantum creation theory 
[16j or a quintessence scalar field cosmology Q. As we 
will discuss later on, e has to be very small (i.e., e <C 1), 
otherwise the radiation density would become much too 
big, contrary to present values. As an order of magnitude 
estimate we expect e ~ p^_Pj.^ fQj- which indeed e <C 1, 
since p-y + p-y <^ p,„. 

Adopting as mentioned above the relation px = —px 
and defining px = A(t)/(87rG), from Eq. (|2.7p we obtain 
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(2.8) 



We assume a power law model for the A{t) function, 
A{t) = B{R/Rq)~'^ , or equivalently in terms of redshift, 
A(z) = B{\ + z)™, where S is a constant. The value of 
this constant is i? = 2>Hq{1 — ^mo), which can be found 
using Eq. (|2.ip at the present epoch, assuming that to- 
day p^ <^ PrrnPx- Thus, the density evolution for the x 
component is given by 



Px/Pcrit = (1 - f^mo)(l + zY 



(2.9) 



where we have defined as Pcnt = f^fQ the present critical 
density of the Universe. If m is positive, then the DE 
slowly decreases as a function of the cosmic time, whereas 
if m is negative the inverse process happens. 

From Eqs. (|2.5p and (|2.6p we derive the evolution laws 
for matter and radiation. 



Pm/Pcrit = ^mo(l + zf 
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where ^-caa ^-nd fi-yo a-re the matter and radiation energy 
densities at z = 0, respectively. 

Since we are interested in the evohition of the radiation 
temperature, it is useful to discuss the extreme case when 
only photons enter in the process (i.e. e = 1). Then 
for the matter density in Eq. (|2.10p the usual evolution 
a (1 + 2)'^ holds, while the radiation besides the usual 
term r2^o(l + z)^''' has also a perturbative term depending 
on 771. Since today ^ Jl-yo ~ 5 x 10~'^ it turns out that m 
has to be extremely small <^ 10~^. Therefore, unless m 
is extremely small or vanishing, DE has to decay mainly 
in matter with possibly some photons as well. Thus the 
condition e ^ 1 has to hold. 



interesting quantity is the deceleration parameter, which 
can be written as 



RR 

q(z] = : — 



(2.16) 



(1 + z)3(m - Sn^o) + 3(m - 2)(1 + z)"'{n^o - 1) 
2(1 + z)3(m - Sano) + 6(1 + z)™(a„o - 1) 

Imposing that q{z) ~ 0, we can determine the transition 
redshift, i.e. the redshift at which the Universe changed 
from a deceleration to an acceleration phase, which is 
given by 
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From this result we see that the larger m is, the earlier 
the Universe changes from deceleration to acceleration. 



C. Hubble and deceleration parameter 

Due to the very small value of Hjq it follows that the 
evolution of the Universe is essentially driven by the DE 
and DM, therefore, from Eq. (|2.ip the following law for 
the Hubble parameter holds 
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which is obviously the same expression found in Ma [20| . 
Recasting Eq. (|2.7p . it is possible to write 

p, + 3H{p., + p, + ^)^0, (2.13) 

which shows that the term Cx contributes to an effective 
pressure 



Pcff = Px 
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(2.14) 



Therefore, we get an equivalent effective DE equation of 
state WaB [201 



Woff 



Px 3 



(2.15) 



If m > then we have Wcs > —1, i.e. our model is 
quintessence-like [3, |^, [l3| , while we have a phantom-like 
[lO| model when m is negative and w^ff < — 1. Another 



D. Thermodynamical aspects and CMB 
temperature evolution 

We follow here the approach outlined in Lima jl5| . 
where he defines a current as iV" := nu" with n being the 
particle number density of the photons or of the DM par- 
ticles. Indeed, there is a current for each of these compo- 
nents. Due to the decaying vacuum the current satisfies 
the following balance equation (one for each component) 



3niH = ipi 



(2.18) 



with i = 7 or DM (7 for the photons and DM for the 
dark matter) and ipi is the corresponding particle source. 
For decaying vacuum models ip-y -f V'da/ is positive and 
related to the rate of change of p^- We can also define 
an entropy current of the form 
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where cr, is the specific entropy per particle (photons, 
DM and in principle also DE). If the DE p^ is constant 
the above entropy current is conserved. The existence of 
a non equilibrium decay process of the vacuum implies 
•Sq > 0, thus an increase of the entropy as a consequence 
of the second law of thermodynamics. In principle the 
second law should be applied to the system as a whole, 
thus including the vacuum component [25|. Assuming 
that the vacuum is like a condensate with zero chemical 
potential pvac it follows from Eulcr's relation 
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(2.20) 



^ The present radiation density is the only cosmological parameter 
accurately measured. The radiation density is dominated by the 
energy in the cosmic microwave background (CMB), and the 
COBE satellite FIRAS experiment determined its temperature 
tobeT = 2.725 ± 0.001 i^ |2l, corresponding to H^o ~ SxlO"^. 



provided Wx = —1, that a^ac = and thus its contribu- 
tion to the entropy current vanishes. Given our assump- 
tions the vacuum plays the role of a condensate carrying 
no entropy. 
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FIG. 1: Radiation temperature in Eq. (|2.33p in terms of m, fimo and 7. Red, blue and green lines are for z — 0.5, z = 1 and 
z = 2, respectively. Left. T as a function of m. 7 = 4/3 and fimo = 0.1 (long-dashed), = 0.3 (continue) and = 0.5 (short- 
dashed). Middle. T as a function of flmo- 7 = 4/3 and m = —0.2 (long-dashed), = (continue) and — 0.2 (short-dashed). 
Right. T as a function of 7. QmO = 0.273 and m = —0.2 (long-dashed), = (continue) and = 0.2 (short-dashed). 



For instance, for quintessence models in the limit where 
the scalar field does not depend on time, and thus its time 
derivative vanishes, one gets w^; = — 1. In the later stages 
of the Universe the time dependence is possibly very weak 
so that Wx = —1 holds up to small corrections. 

The equation for the particle number density of radia- 
tion component is given by Eq. (|2.18p with i = 7. Using 
Gibbs law and well-known thermodynamic identities, fol- 
lowing the derivation given in the paper by Lima et al. 
[l7| . one gets (see also [ll.[l5|) 
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(2.21) 

To get a black-body spectrum the second term in 
brackets in Eq. p.2ip has to vanish, thus 
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Now, following the previous discussion on e and aiming 
to be very general, we set e = ti—Ej-l^ where e is a new 
parameter and insert it into Eq. (j2.26p . Taking the sum 
of Eqs. (HH) and (g^ we find 



8^G{ptot+Ptot) ^ SnGpm = 2I2 -2| = -2H . (2.27) 



Finally, we obtain the expression 
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which we can integrate 
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Thus, Eq. (P?^ becomes 



With 
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(7—1) one obtains 
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Using the equation for the particle number conserva- 
tion Eq. ([2l8| into Eq. ([2:24| leads to 
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%dt = (7 - 1) 
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(2.29) 



where to denotes the present time and <i some far instant 
in the past. Indeed, if A vanishes and 7 = 4/3 one gets 
the usual dependence T{t) = ^ ^,^!^ ''' for a radiation 
fluid. To carry out the integration of the first term on 
the right hand side it is useful to perform a change of 
variable from t to z and accordingly ^ = jyTITT)' This 
way we get (with 21 corresponding to the time ti and 
zq = corresponding to to present time) 



, Tiz = 0) ^, ^,, i?(z = 0) 
In '^, ' + 37 - l)ln p / 
T(zi) R{zi) 

(7-I) n A^f~ 

2 7o H'H{1 + z) 



(2.30) 



dz 



With Eqs. (11221) and ^^TI^i we get 



where ' denotes derivative with respect to z. 



As next we insert H{z) and its derivative as taken from 
Eq. (PJ^ into Eq. ([^30)1 and integrate it, to get (setting 

zi = z) 



T{z) = To 
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where 



A = ln{{m - Sano) + m(l + z)™-^^^^^^ _ ^^^ (2.32) 

-ln{{m- 3)fi,„o)] ■ 



We can also write Eq. p.3ip as 



T{z) = ro(i + z) 



3(7-1) 



(m ^ 3r!.„o) + m(l + z)"-^(r!.„o - 1) 
(m - 3)f2,„o 



(2.33) 

?(7-l) 



Wc inserted in the exponent of Eq. (|2.3ip the cxphcit 
form of B, thus getting as exponent in the above Eq. 
e(7 — 1). Hereafter, we wiU set e = 1. Clearly e and m 
are not independent, we checked using the temperature 
redshift data that if e is bigger (~ 10 or more), then m has 
to be extremely small consistently with what mentioned 
in section II. B (as it would lead to a too high production 
of photons in the DE decay). On the other hand, if e 
gets smaller (e.g., ~ 0.1) m gets bigger (^ 0.2) and ac- 
cordingly Weff, moreover m would be poorly constrained, 
since the uncertainties would then be very high. But from 
the other data (without the the temperature ones) there 
are already stringent limits on m and thus this way one 
could get lower limits on e, under the assumption that 
DE decays also in photons. 

Notice that for z = we have T(0) ~ To, whereas for 
TO = the expression in the parenthesis is equal to 1 
and thus T{z) ~ ro(l + z)'^^'^~'^\ which for the canonical 
value of 7 = 4/3 reduces to the standard expression. 

In Fig. [TJ the temperature relation from Eq. (|2.33p 
is discussed in terms of the model parameters and red- 
shift. At fixed redshift, the temperature is a decreas- 
ing function of m, i.e. larger m means a colder CMB 
radiation temperature and when compared with the un- 
perturbed case with m = 0, then a decaying/increasing 
model predicts colder/hotter temperatures. The temper- 
ature is sensitive to Clyno mainly at r^mo ^ 0.1 — 0.2 and 
less at larger flmOj while it is a strong function of 7, in- 
creasingly as a function of redshift. 



III. FITTING PROCEDURE AND DATA 

To constrain the parameters of the model we use a set 
of different kinds of measurements. Our data sets include 
the measurements of CMB temperatures from high red- 
shift quasars and SZ effect, high quality "Unionll" SN 
la data, baryon acoustic oscillation measurement from 
the Sloan Digital Sky Survey, the shift parameter from 



WMAP three years results and 9 observational H{z) 
data. To break the degeneracies between the parameters 
and explore the power and differences of the constraints 
for these data sets, we use them in several combinations 
to perform our fitting. 

In order to constrain the model parameters we maxi- 
mize the likelihood function £(p) ~ exp [— jX^] 1 where p 
denotes the set of model parameters and x^ is a suitable 
merit function'^. The isolikelihood (or the isox^) con- 
tours provide constraints on the parameter space. The 
68% confidence levels (CL) are obtained by imposing 



Ax' 



A Ari 



1 and 2.3 for 



1, 2 free param- 



eters, where Xmin is the minimum of x^ function. The 
90% CL, is given by Ax^ = 2.71 and 4.61 for Up = 1, 
2. Finally, the 95% CL is given by Ax^ = 4 and 6.17 
for Hp — I, 2. In order to give a quantification of the 
errors on a given parameter we can follow different ap- 
proaches. When rip = 1, then the error on the parameter 
is determined simply adopting the above conditions on 
Ax^. On the contrary, if ?ip = 2 (which is one of the 
cases we will investigate in the paper), to constrain a 
given parameter pi , we rely on the marginalized function 
defined as £pj (pi) oc J dp2C{p), which is normalized to 
1 at the maximum. The value of pi corresponding to the 
maximum of such a function is chosen as our best fitted 
value'^ and the CLs are determined by applying to this 
marginalized likelihood the conditions above on Ax^ for 
Up = 1. 

In the course of the paper, for sake of simplicity we 
will report only the 68% CL for the listed best fitted 
parameters. 



A. Temperature measurements (T) 

To test the temperature evolution for the radiation 
component, we rely on the CMB temperatures derived 
from the absorption lines of high redshift systems and 
the ones from SZ effect in clusters of galaxies (wc will 
collectively quote as Tcmb, hereafter). At high redshift 
the CMB temperature is recovered from the excitation of 
interstellar atomic or molecular species that have tran- 
sition energies in the sub-millimetre range and can be 
excited by CMB photons. When the relative population 
of the different energy levels are in radiative equilibrium 
with the CMB radiation, the excitation temperature of 
the species equals that of the black-body radiation at 
that redshift, providing one of the best tools for deter- 
mining the black-body temperature of the CMB in the 
distant Universe [27H34| . In Jetzer et al. [l|, we adopted 
a sample of 5 QSO adsorption measurements, which is 
now updated to 9 after the recent measurements reported 



^ This is equivalent to a x^ minimization. 

^ Note that for asymmetric CLs, the x^ minimum and the maxi- 
mum of the marginalized likelihood can be different. 



in Notcrdaeme et al. 12J]. In summary we have 4 data 
points from the analysis of the fine structure of atomic 
carbon (AC) and 5 measurements based on the rotational 
excitation of CO molecules (CO) |24| . 

At lower redshift we use the measurements from the 
SZ effect. During passage through a cluster of galaxies 
some of the photons of the CMB radiation are scattered 
by electrons in the hot intracluster medium. This im- 
print was first described by SZ [35|. Thus, spectral mea- 
surements of galaxy clusters at different frequency bands 
yield independent intensity ratios for each cluster. The 
combinations of these measured ratios permit to extract 
the cosmic microwave background radiation (see Fabbri 
et al. 1361') . We will rely on the data compilation in Luzzi 
et al. j23l |. which have analyzed the results of multifre- 
quency SZ measurements toward several clusters from 5 
telescopes (BIMA, OVRO, SUZI II, SCUBA and MITO). 

We will match the observed Tcmb with the theoretical 
expression Tth, which we have derived in Eq. (j2.33p . by 
minimizing the following merit function 
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(3.1) 



where <JcMB,i is the error on the temperature estimates 
and Ntcmb ~ 22 is the number of available observa- 
tional data. 



B. High quality Supernovae la data set (SN) 



where 
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(3.5) 
(3.6) 
(3.7) 



This new function does not depend on /^Oi allowing us to 
drop the contribution from the Hubble constant. 



C. Baryon Acoustic Oscillation (A) 

In the large-scale clustering of galaxies, the baryon 
acoustic oscillation signatures could be seen as a stan- 
dard ruler providing the other important way to con- 
strain the expansion history of the Universe. We use the 
measurement of the B AO peak from a spectroscopic sam- 
ple of 46,748 luminous red gala xies (LRGs) observations 
of SDSS to test cosmology [421, which gives the value of 
A = 0.469(nf/0.98)-°-35 ± 0.017 at zbao = 0.35 where 
Us — 0.96 [43[. The expression of A can be written as 
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(3.8) 



The most important candle we use is the type la su- 
pernovae (SN). We adopt the Unionll dataset discussed 
in AmanuUah et al. [33], which consists of Nsn = 557 
datapoints from z = to z = 1.4, compiled after the com- 
bination of different datasets and the consequent appli- 
cation of various selection cuts to create a homogeneous 
and high signal-to-noisc sample. 

The data points for SN are given in terms of distance 
modulus Hobs = m—M, where m and M are the apparent 
and absolute magnitude, respectively. The theoretical 
distance modulus is given by 



fithiz) == 5 logio Dl{z) + Ho 



(3.2) 



where ^o = 42.38 — 5 logj^g h and Dl(z) is the luminosity 
distance at the redshift z. The x^ function to be mini- 
mized is 
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where asN.i is the error on Hobs,i- The parameter ^o 
is a nuisance parameter which depends on the Hubble 
constant. One can perform a standard marginalization 
on fiQ- Otherwise, following |38l - l4l| . it is easy to check 
that the x^ in Eq. p.3p is equivalent to the following 
function 



~2 _ A ^' 
XsN — ^ pi") 



(3.4) 



which is evidently independent on Ho, and the relative 
X^ function is 



2 _ /^-0.469(ns/0.98)-°-35N^2 
Xbao - [ 5;^Y7 



D. CMB Data: the shift parameter (R) 



(3.9) 



The measurement of CMB anisotropics represents a 
powerful tool to constrain cosmological parameters, how- 
ever, using the full data of CMB is time consuming, thus 
as an alternative it is common to rely on the measure- 
ment of the shift parameter R. The CMB shift parameter 
may provide an effective way to constrain the parameters 
of DE models since it has the very large redshift distribu- 
tion, which allows to constrain the evolution of DE very 
well. The shift parameter R which is derived from the 
CMB data takes the form as 



R 



Vn 



mO 



^CMB 



dz' 



H{z')/Ho' 



(3.10) 



where is zcmb = 1090 and the observed value for Eq. 
(|3T0)) has been updated to i? = 1.71 ± 0.019 from 
WMAP5 0. The x^ function is 



XcMB 



i?-1.71 
0.019 



(3.11) 



IV. RESULTS 
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FIG. 2: Cosmic microwave background temperature as a func- 
tion of the redshift. The black points with bars are the fuU 
collection of measurements from Luzzi et al. [23| and Noter- 
daeme et al. [2J|. The gray arrows represent the upper limits 
derived from the analysis of atomic carbon (see [2J| for de- 
tails). The black line is the best fit result (m = 0.03), while 
the gray region is the la uncertainty. The red line is the best 
fit recovered from Ma [20|]. The inset panel show a magnified 
vision of the higher redshift region of the plot. 



E. Observational H{z) Data (OHD) 

By using the differential ages of passively evolving 
galaxies determined from the Gemini Deep Deep Survey 
(GDDS) and archival data [4J|, Simon et al. determined 
H(z) in the range < 2 < 1.8 [4^. The 9 observational 
Hobs,i datapoints can be obtained from [45l - l47l |. The x^ 
statistics for these H{z) data is 



XoHD ~ 7 , 



i=l 



2 

^OHDA 



where aoHD,i is the error on logj^Q Hobs,i- Following the 
same procedure adopted in Scc lIIIBl to marginalize with 
respect to Hq, we replace the x^ in Eq. p.l2p with Eq. 
(lOl). where 



^-E 



1=1 
9 



(logio H{z^,Ho = 1) - logio Hob 

2 
^OHDA 



, (3.13) 



B^Y^ (logio H{z,,Ho = 1) - logio Hobs,i) ,^ ^^. 



4=1 

9 



2 






We discuss in this section the best fitted values for 
our model parameters. We first concentrate on the 
temperaturc-redshift relation, deriving constraints on 
both 771 and 7. In order to constrain the present mat- 
ter density J7nio and m we add the set of observational 
probes we have listed in the previous section, which are 
primarily linked to the Hubble parameter H{z), and thus 
much more sensitive to matter component. 



A. Constraints from T{z) — z relation 

As a primary test, following the same line in Jetzer 
et al. [l|, we have compared the CMB temperature pre- 
dicted (see Eq. (|2.33p ). with the updated collection of 
multi-redshift measurements of Tcmb we have discussed 
in the previous section. We set Tq = 2.725 K, which 
is quite well determined in the literature [48|, and the 
matter densit y f^ mo = 0.273 to the value inferred in Ko- 
matsu et al. [4j]. As we have shown in Sect. IIIDl the 
temperature is not sensitive to changes in J7niO in the 
region where J7mO >0-l ~ 0.2 in which it is constrained 
to lie from other probes. Thus, we are sure that our 
estimates are robust and in this section we will not dis- 
cuss the constraints on f7,„o coming from the fitting of 
temperaturc-redshift relation. If we take 7 = 4/3, then 
we find m = 0.03^o'Qg, which is lower than the estimated 
value of m = 0.09 ±0.10 in [l|, but fully consistent within 
uncertainties, and also pretty consistent with to = 0. 
In Fig. [2] the temperature measurements (together with 
some upper limits) are shown, and our best fitted result 
is plotted and compared with the to, = —0.09 result in Ma 
[201 . The value we have found corresponds to an effec- 
tive equation of state Weff = —0.99 ±0.03, consistent with 
Weff = — 1, and the transition redshift is zt = 0.78 ±0.08. 
In order to check the impact of redshift distribution we 
separate the data in two redshift bins with z < 0.6 (SZ 
data only) and z > 0.6 (QSO adsorption lines only), find- 



(3.12) ing the best fitted values tti = 0.12 



+°-i2 and 

0.13 ^^^ 



-0.05 



+0.12 

-0.14J 



respectively. Although the uncertainties are very high 
and no statistically relevant conclusion can be reached, 
these results give some indications of a mild trend with 
lower redshift data preferring a DE decaying into matter 
and radiation, while data at z > 0.6 point to an op- 
posite behavior. These results could be interpreted in 
a different way, in fact, the differences found could be 
due not to the different redshift coverage, but to some 
particular biases in the two kind of observations, SZ vs 
QSO adsorption lines. An indication of this suggestion 
comes if we divide the sample in three subsamples: 1) SZ 
data, 2) the data from the analysis of the fine structure 
of atomic carbon (AC), and 3) the measurements based 
on the rotational excitation of CO molecules in [2J|. If 
we fit the model to the combined SZ±AC and SZ±CO 



(3.15) samples we find 



0.11 



-fO.ll 

0.11 



and TO = 0.03 



-1-0.09 
-0.10' 



re- 



spectively. Because of the larger measurement errors. 
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FIG. 3: The 68%, 90% and 95% confidence limit contours in 
the m— 7 plane from the fitting of temperature measurements. 
The red, blue and green contours are relative to the fit using 
all the sample, data for z < 0.6 and the ones for z > 0.6, 
respectively. The points are the best fitted values, while boxes 
are the results when 7 = 4/3. 



the AC data affect very little the results when only SZ 
measurements are adopted, and the result for the seeond 
sample shows that SZ and CO data mainly constrain m, 
and AC simply gives a tiny reduction on the errors. Of 
course, larger data samples would be needed to provide 
definitive answers. 

Adopting a constant value for the ratio ip^/in^H — f3 
Lima et al. [131 have found the simple relation. 



Tiz) = To{l + zy-P 



(4.1) 



If we fit ah the sample we find ^ = -0.002 ± 0.03, 
while for the low and high redshift subsamples we have 
/? = 0.06 ± 0.08 and -0.01 ± 0.03, respectively. These 
results are qualitatively consistent with what found in 
Noterdaeme et al. [2^, and points to a similar trend as 
the one discussed above^. 

If 7 is left free to change, we obtain the contours shown 
in Fig. [3] In all the redshift samples, for the best fit- 
ted value it turns out that 7 > 4/3 and m is systemati- 
cally more positive. Adopting the whole sample, the CL 
contours are broad, and from the marginalization with 
respect to the second parameter we find 7 ~ 1.35 



and 



0.25 



+0.23 
-0.17' 



+0.03 
0.03 



which corresponds to an effective 



equation of state Wcs = —0.92 ± 0.07 and the transition 
redshift is zt ~ 1.1 ± 0.6. When the two subsamples 



* Although we use the same datasample as in |24| . we find some 
difi^erences (although very minor), which could be possibly due 
to the way the asymmetric errors in QSO absorption line data 
are accounted for. In particular, we have adopted as error in the 
fit the average of the two errors. 



TABLE I: Maximum likelihood parameter and la uncertain- 
ties of 771 for QmO = 0.273 and for m and f^mo when this last is 
left free to vary. The legend of the symbols is: T = Temper- 
ature, SN = Supernovae la, A — Baryon Acoustic oscillation 
parameter, R — Shift parameter, OHD — Observational H{z) 
data. 



Model 


n„bO = 0.273 


f^mO 


free 




m 


f^mO 


m 


SN 


0.01+^1? 


0.251°:^^ 


-0.2^^7 


SN-FT 


0-03tVol 


0.281°:°^ 


0.03l°;J° 


SN-hA 


0.01+Si? 


90+0.02 


0.03l«:^| 


SN-^A-fR 


-OMt°f, 


0.27t°:S? 


-0.041H1 


SN-l-A-fR+T 


-o.oslH^ 


0.27t°:°? 


-o.oslH^ 


SN-FA-hOHD 


-omtl;\t 


r, 97+0.03 

"•^'-0.02 


-o.o2l«:?^ 


SN-hA-hOHD+R 


~OM+_°f, 


r. r,y+0.02 

'J-^i -0.01 


-OM+Jof, 



arc adopted, wide confidence contours are found, partic- 
ularly for the z > 0.6 sample, for which the contours at 
very low m are not closed. We obtain 7 = 1.3lg'^ and 

1.26to;^i, while m = O.StaJ and Q.6t%, respectively for 
the low and high-z samples^. 



B. Constraints from independent measurements 

If we set 17,110 = 0.273, we can give some further con- 
straints (shown in Table H]) on m, using the other ob- 
servational probes listed in Sect. IIIII together with tem- 
perature measurements. We note that SNs alone are not 
able to constrain the value of m, producing very high un- 
certainties, also larger than the ones obtained using the 
temperature measurements alone. The Baryon Acous- 
tic Oscillation parameter does not help, while tempera- 
ture measurements will allow to reduce the uncertainties. 
These probes produce slightly positive values for to, while 
adding the Shift parameter and/or the H{z) data give 
negative to ~ —0.03, —0.04. Anyway the uncertainties 
remain very large, with the only exception of the results 
when the Shift parameter is used in the fitting, m being 
in this case greatly constrained to m 
Table H] for details. 

When fimo is left free to vary, we find the best fitted 
values and CL contours shown in Table |T] and Fig. H] 
Remarkably, the best fitted values of f^nio are perfectly 
in a gree ment with the values found in independent works 
[3 71 . |43| . Consequently, it is not surprising that the re- 
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We notice that due to the particular form of the CL contours the 
best fit quantities derived from the maximum of the marginalized 
likelihood can be different from the x^ minimum. In fact, we 
find 7 = 1.6 and 1.4, while m = 0.6 and m = 0.4 for the two 
subsamples, which differ from the ML values reported in the text, 
but consistent within errors, due to the large uncertainties. 
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FIG. 4: The 68%, 90% and 95% confidence limit contours in the fimo —rn plane. See the legends for the meaning of the symbols. 



covered values for ttt,, except for the case when only SNs 
arc fitted, are still quite consistent with the estimates 
obtained when riuio is fixed (sec Table H]). If we con- 
sider the fit with SN+A+OHD+R, then our best fitted 
parameters correspond to an effective equation of state 
Woff = —1.01 ± 0.01, in the phantom regime but fully 
consistent with WcS = — 1, and the transition redshift is 
ZT = 0.72 ± 0.04. 

Adopting a standard DE model with an equation of 
state Px = WxPx, and Wx free to vary, we have performed 
the fit, finding that Wx ^ —1, consistently with our re- 
sults and still pointing to a cosmological constant as the 
best description for DE. 



V. CONCLUSIONS 

We have presented the properties of a variable DE 
model, following the approach in Jetzer ct al. [l| (see 
also [ISl - llSl . [20] )■ The model relies on the assumption 
that the DE component is not conserved and exchanges 
particles with both DM and radiation components. In 
particular, we show that in order to have a viable model, 
DE needs to decay besides into radiation also in matter. 
Motivated by these considerations, we have concentrated 
the analysis on the case when only a small fraction is 
exchanged with radiation. We have dedicated partic- 
ular attention to the thermodynamic properties of the 
model, discussing the theoretical relation between the ra- 
diation temperature and redshift and we have matched 
it to the most updated collection of CMB temperature 
measurements [l|, [2j, [2J] , consisting of both SZ data and 
high-redshift QSO adsorption line observations. First, we 
have constrained the model by setting 7 = 4/3 and find- 
ing m = O.OSJ^QQg, consistent with our previous estimate 



within uncertainties, still consistent with the standard 
case of a T ex (1 + z), i.e. with an effective equation of 
state Wcff = —1. When 7 is left free, then m = 0.25J^Q'^y. 
Although the large contours indicate that is not possi- 
ble to find statistically relevant departures from a stan- 
dard cosmological constant, we find that a decaying DE, 
with an effective equation of state Wcs > — 1 is preferred. 
Then, for the first time, we test this kind of model com- 
bining both CMB temperature measurements and dif- 
ferent observational data sets, like Supernovae la, CMB 
and large-scale structure data, etc. We find that il^o is 
almost independent on the combination of data we use. 



and in the best case is constrained to i}„ 



0.27 



+0.02 

-0.01- 



Moreover, while temperature measurements and SN data 
tend to furnish m >0, consistently with a decaying DE 
model, other datasets, like CMB data prefer the oppo- 
site situation with to ^ 0, i.e. a phantom-like model with 

Wcff < -1- 

Although the present data do not allow to find 
strong discrepancies with a classical cosmological con- 
stant model, we think that future surveys at high red- 
shift could collect further measurements for CMB tem- 
perature, which could help us to further constrain the 
temperature-redshift relation. This kind of information, 
together with larger and higher quality samples of Su- 
pernovae and better CMB data or other standard candles 
like Gamma ray burst (GRBs) can allow to give indepen- 
dent constraints on both matter density i7,„o and effective 
equation of state Wcs- 

We notice that if DE docs not decay into radiation 
(corresponding to e = e = and thus to, is no longer 
constrained) then the CMB temperature will scale in the 
standard way. Clearly, this would imply that if DE de- 
cays, this has to be into DM only. On the other hand 
a deviation of the CMB temperature from the standard 
scaling could be interpreted as DE decaying also into ra- 



10 



diation. In which case with Eq. (|2.33p one can determine 
m and/or i and thus get some insights on the decay mode 
of DE into radiation. Future data on the CMB temper- 
ature will allow to shed light on this important issue. 
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